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The large-scale organization of the world economies is exhibiting increasingly levels of local het- 
erogeneity and global interdependency. Understanding the relation between local and global features 
calls for analytical tools able to uncover the global emerging organization of the international trade 
network. Here we analyze the world network of bilateral trade imbalances and characterize its over- 
all flux organization, unraveling local and global high-flux pathways that define the backbone of the 
trade system. We develop a general procedure capable to progressively filter out in a consistent and 
quantitative way the dominant trade channels. This procedure is completely general and can be 
applied to any weighted network to detect the underlying structure of transport flows. The trade 
fluxes properties of the world trade web determines a ranking of trade partnerships that highlights 
global interdependencies, providing information not accessible by simple local analysis. The present 
work provides new quantitative tools for a dynamical approach to the propagation of economic 
crises. 



I. INTRODUCTION 

The term "globalization" , when applied to the inter- 
national economic order, refers to the presence of an in- 
tricate network of economic partnership among an in- 
creasing number of countries [T]. In this context, the 
International trade system, describing the fundamental 
exchange of goods and services, plays a central role as 
one of the most important interaction channels between 
states P]. For instance, it broadly defines the substrate 
for the spreading of major economic crises [3J 4, 5 , such 
as the 1997 Asiatic crisis [SJ EJ |7] which shows how eco- 
nomic perturbations originated in a single country can 
somehow propagate globally in the world. Moreover, 
commercial trade flows are indeed highly correlated with 
other types of cross-country economic interactions (flows 
of services, financial assets, workers) and so stand as a 
good indicator for more general economic relations [5]. 
The International trade system as an independent extract 
of the world economy is therefore still a partial view of the 
whole system; a complete description would consider the 
feedback mechanisms that operate between international 
trade imbalances and other economic variables such as 
investment, debt, or currency prices. On the other hand, 
the study of the International trade network in a system's 
perspective represents a necessary first step before pro- 
ceeding with a subsequent more integrative investigation 
and has proven to be successful in providing insight into 
some of its global properties. 

The large size and the entangled connectivity pattern 
characterizing the international trade organization point 
out to a complex system whose properties depend on its 
global structure. In this perspective, it appears natu- 
ral to analyze the world trade system at a global level, 



every country being important regardless of its size or 
wealth and fully considering all the trade relationships. 
A convenient framework for the analysis of complex in- 
terconnected systems is network analysis jH [10] . Within 
this approach, countries are represented as nodes and 
trade relationships among them as links. Such visual- 
izations of bilateral trade relations have been used in 
recent years to help analyze gravity models jTTJ [TJ] , of- 
ten proposed to account for the world trade patterns and 
their evolution 13J. While the first attempts to study the 
trade system as a complex network have successfully re- 
vealed a hierarchical organization [TH [T5J [TB] , these stud- 
ies focused on topological aspects neglecting fundamental 
components, such as the heterogeneity in the magnitude 
of the different bilateral trade relations and their asym- 
metry. These are essential issues in the understanding of 
the interplay between the underlying structure and the 
principles that rule the functional organization of the sys- 
tem. 

Here we tackle the quantitative study of the world 
trade network by implementing the trade flux analysis 
at a global scale. To this end, we construct the weighted 
directed network of merchandize trade imbalances be- 
tween world countries. In this representation, each coun- 
try appears as a node and a directed link is drawn among 
any pair whenever a bilateral trade imbalance exists, 
i.e., whenever bilateral imports does not balance exports. 
The direction of the arrow follows that of the net flow of 
money and it is weighted according to the magnitude of 
the imbalance between the two countries. More precisely, 
we define the elements Eij that measure the exports of 
country i to country j and the elements Iij that measure 
the imports of country i from country j. The trade im- 
balance matrix is therefore defined as Tj ? - = En — In and 
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FIG. 1: Measuring local inhomogeneity in fluxes, a, concep- 
tual representation of the link construction process, b and c, 
local inhomogeneity for incoming (b) and outgoing (c) con- 
nections measured by kY(k) as compared to the null model. 
The diagonal line corresponds to the maximum possible in- 
homogeneity, with only one connection carrying all the flux. 
The line kY(k) — 1 is the maximum homogeneity, with all 
the fluxes equally distributed among the connections. The 
area depicted in grey corresponds to the average of kY(k) un- 
der the null model plus two standard deviations. The solid 
lines are the best fit estimates which e ive kinY (^kinj ~ kin 
and k utY(k 0U t) ~ k^t- The inset in (c) sketches a pathway 
through a country arising from strong local inhomogeneity in 
incoming and outgoing connections. 



measures the net money flow from country j to country 
i due to trade exchanges. Since E^ — Iji and = Eji, 
T is an antisymmetric matrix with = —Tji, and a di- 
rected network can be easily constructed by assuming a 
directed edge pointing to the country with positive bal- 
ance. The network of the net trade flows is therefore 
defined in terms of a weighted adjacency matrix F with 
Fij =| Tij | = | | for all i,j with T i} < 0, and F %] = 
for all i,j with > (see Fig. la for a pictorial de- 
scription) . 

By using the above procedure we constructed the net- 
work of trade imbalances by using the data set which 
reports the annual merchandize trade activity between 
independent states in the world during the period 1948- 
2000, together with the annual values of their Gross Do- 



mestic Product per capita and population figures (1950- 
2000) [17], [25]. The time span of the data set allows us 
to study the change of trade flow networks with yearly 
snapshots characterizing the time evolution of the trade 
system. The most basic topological characterization of 
each country within the network is given by the number 
of incoming and outgoing links, k in and k out respectively, 
which inform us about the number of neighboring coun- 
tries that contribute positively and negatively to the net 
trade imbalance of the country under consideration. A 
precise assessment of the country trade balance cannot 
however neglect the magnitude of the fluxes carried by 
each trade relation. This information can be retrieved 
summing up all the weights of the incoming or outgo- 
ing links, which give us the total flux of money due to 
trade entering to or leaving from the country of inter- 
est. In the network literature, these two variables are 
called incoming and outgoing strength and are denoted 



by 



J2i F ij and s °j ut = Si-Pjii respectively 



The total trade imbalance of a country can then be com- 
puted as Asj = s" 1 — s° ut . Depending on As., , countries 
can be then defined as net consumers and net produc- 
ers. Net producers export more than they import, the 
total outcome being a trade surplus which corresponds 
to Asj > 0, whereas net consumers export less than they 
import, the total outcome being a trade deficit which is 
indicated by Asj < 0. Since one incoming link for a given 
country is always an outgoing link for another, the sum 
of all the countries' trade imbalances in the network must 
be zero. While the local balance is not conserved, we are 
therefore dealing with a closed system which is globally 
balanced (the total flux is conserved). Merchandizes, or 
equivalently money, flows in the system from country to 
country with the peculiarity that there is a global flow of 
money from consumer countries to producer ones. 



II. LOCAL HETEROGENEITY AND 
BACKBONE EXTRACTION 

The obtained networks show a high density of connec- 
tions and heterogeneity of the respective fluxes among 
countries. Indeed, as the number of countries increases, 
so does the average number of trade partners, as well as 
the total flux of the system, which is seen to grow pro- 
portional to the aggregated world Gross Domestic Prod- 
uct [TJ] . The overall flux organization at the global scale 
can be characterized by the study of the flux distribution. 
A first indicator of the system heterogeneity is provided 
by the probability distribution P{Fij) denoting the prob- 
ability that any given link is carrying a flux Fij. The 
observed distribution is heavy-tailed and spans approx- 
imately four orders of magnitude [TJ]. Such a feature 
implies that only a small percentage of all the connec- 
tions in the network carry most of its total flow F and 
that there is no characteristic flux in the system, with 
most of the fluxes below the average and some of them 
with a much higher value. This is however not totally 
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TABLE I: Sizes of the backbones. Percentage of the original 
total weight F, number of nodes N and links E in the 1960 
and 2000 imbalance networks that remain in the backbone as 
a function of the significance level a. 
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88 


100 25 
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98 25 
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83 
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79 
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69 


92 9 


75 
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unexpected since a large scale heterogeneity is a typical 
feature of large-scale networks. In addition, the global 
heterogeneity could just be due to differences in the sizes 
of the countries, in their population and in their respec- 
tive Gross Domestic Product. More interesting is there- 
fore the characterization of the local heterogeneity; i.e. 
given all the connections associated to each given coun- 
try, how is the flux distribution for each of them. 

A local heterogeneity implies that only a few links 
carry the biggest proportion of the country's total in-flow 
or out-flow. Interestingly, such a heterogeneity would de- 
fine specific pathways within the network that accumu- 
late most of the total flux. In order to asses the effect 
of inhomogcncitics at the local level, for each country 
i with fc incoming or outgoing trade partners we calcu- 
late na e] 

k 

kY i (k) = k'£p 2 ij , (1) 
j=i 

where k can be either k in or k out in order to discern be- 
tween inhomogeneities in incoming and outgoing fluxes, 
and where the normalized fluxes of node i with its neigh- 
bors are calculated as pij = Fji/sl" 1 for incoming con- 
nections and as p^ — Fij/s° ut for the outgoing ones. 
The function Yi(k) is extensively used in economics as a 
standard indicator of market concentration, referred as 
the Herfindahl-Hirschman Index or HHI [55J H3] > and it 
was also introduced in the complex networks literature 
as the disparity measure [23]. In all cases, Yj(fe) char- 
acterizes the level of local heterogeneity. If all fluxes 
emanating from or arriving to a certain country are of 
the same magnitude, kYi(k) scales as 1 independently 
of fc, whereas this quantity depends linearly on k if the 
local flux is heterogeneously organized with a few main 
directions. Increasing deviations from the constant be- 
havior are therefore indicating heterogeneous situations 
in which fluxes leaving or entering each country are pro- 
gressively peaked on a small number of links with the 
remaining connections carrying just a small fraction of 
the total trade flow. On the other hand, the deviations 
from the constant behavior have to be expected for low 
values of k and it is important to compare the obtained 
results with the deviations simply produced by statistical 
fluctuations. To this end, we introduce a null model for 



the distribution of flows among a given number of neigh- 
bors in order to assess, in a case by case basis, whether 
the observed inhomogeneity can just be due to fluctua- 
tions or it is really significant. 

The null model with the maximum random homogene- 
ity corresponds to the process of throwing k — 1 points 
in a [0, 1] interval, so that the interval ends up divided in 
k sections of different lengths representing the different 
values assigned to the k variables in the random case. 
It can be analytically proved that the probability that 
one of these variables takes a particular value x depends 
on k and is 

Prob{x < p^ < x + dx} = (k - 1)(1 - x) k ~ 2 dx. (2) 

This probability density function can be used to calculate 
the statistics of kY NM (k) for the null model. Both the 
average and the standard deviation are found to depend 
on k: 

2k 

(kY NM (k)) = k(Y NM (k)) = — — (3) 



so that each node in the network with a certain in or 
out degree should be compared to the corresponding null 
model depending on the appropriate k. 

In Fig. 1, we show the empirical measures along with 
the region defined by the average value of the same quan- 
tity kY{k) plus two standard deviations as given by the 
null model (shadowed area in grey) . For a homogeneously 
random assignment of weights, this quantity converges to 
a constant value for large fc, which is clearly different from 
the observed empirical behavior. Most empirical values 
lie out of the null model domain, which proves that the 
observed heterogeneity is due to a well definite ordering 
principle and not to random fluctuations. 

The direct fit of the data indicates that both in and out 
fluxes follow the scaling law fcYi(fc) oc fc' 3 with /3j n = 0.6 
for the incoming connections and ft ou t =0.5 for the out- 
going ones (see Fig. 1). This scaling represents and in- 
termediate behavior between the two extreme cases of 
perfect homogeneity or heterogeneity but clearly points 
out the existence of strong local inhomogcncitics. The 
emerging picture is therefore consistent with the exis- 
tence of major pathways of trade flux imbalances (thus 
money) that enters the country using its major incom- 
ing links and leaves it through its most inhomogeneous 
outgoing trade channels (see inset in Fig. lc). 

The analysis of the local inhomogeneities in the trade 
fluxes prompts to the presence of high-flux backbones, 
sparse subnetworks of connected trade fluxes carrying 
most of the total flux in the network. This backbone 
is necessarily encoding a wealth of information being the 
dominating structure of the trade system. It is also worth 
remarking that the local heterogeneity is not just encoded 
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FIG. 2: Backbone of the world trade system. Snapshots of the a — 10 -3 backbone of the world trade imbalance web for the 
years 1960 and 2000. Notice that the most central economies are depicted at fixed positions to make both graphs more easily 
comparable. 



in high flux links in terms of their absolute intensities, but 
also takes into account the local heterogeneity by com- 
paring the strength of the fluxes associated to a given 
country with its total strength. It is then interesting to 
filter out this special links and provide snapshots of the 
trade system backbone. This can be achieved by com- 
paring the link fluxes with the null model used for the 
calculation of the disparity in a pure random case. The 
same approach allows us the calculation for each connec- 
tion of a country i of the probability a^- that its normal- 
ized flux value pij is due to chance. Along these lines, we 
can identify highly inhomogcncous fluxes as those which 
satisfy 

fPti 

onj = 1 - (k - 1) / (1 - x) k ~ 2 dx < a, (5) 
Jo 

where a is a fixed significance level. Notice that this ex- 



pression depends on the number of connections of each 
country, k. By choosing a global threshold for all coun- 
tries we obtain a homogeneous criteria that allows us 
to compare inhomogeneitics in countries with different 
number of connections and filter out links that carry 
fluxes which can be considered not compatible with a 
random distribution with an increasing statistical confi- 
dence. The backbone is then obtained by preserving all 
the links which beat the threshold for at least one of the 
two countries at the ends of the link while discounting 
the rest. By changing the significance level we can fil- 
ter out the links focusing on progressively more relevant 
heterogeneities and backbones. 

An important aspect of this new filtering algorithm is 
that it does not belittle small countries and then, it offers 
a systematic procedure to reduce the number of connec- 
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tions without diminishing the number of countries and 
choosing the backbone according to the amount of trade 
flow we intend to characterize. It provides a quantitative 
and consistent way to progressively identify the relevant 
flow backbone once the level of statistical confidence with 
respect to the null case is fixed, or instead the total flow 
we want to represent in the system. Indeed, it is re- 
markable that when looking at the network of the year 
2000 one finds that the a — 0.05 backbone contains only 
15% of the original links yet accounting for 84% of the 
total trade imbalance. Most of the backbones form a 
giant connected component containing most of the coun- 
tries in the network, and only for very high values of 
the confidence level, defining a sort of super-backbones, 
some disconnected components appear and the number 
of countries starts to drop. In this respect, the a = 0.01 
backbone seems to offer the best trade-off since it keeps 
nearly all countries, 75% of the total trade imbalances, 
and one order of magnitude less connections than in the 
original network (see Table 1). 

The backbone reduction is extremely effective in sort- 
ing out the most relevant part of the network and can 
be conveniently used for visualization purposes. For the 
sake of space and reproduction clarity, we report the 
backbones corresponding to a — 10" 3 , still accounting 
for approximately 50% of the total flux of the system. 
Fig. 2 shows two snapshots of such backbones for 1960 
and 2000. These high-flux backbones evidence geographi- 
cal, political and historical relationships among countries 
which affect the observed trade patterns. For instance, 
the trade of US with its geographically closer neighbors 
and also continental neighbors, the case of Russia and the 
former Soviet republics, or France and its former colonics, 
the lack of strong trade relations between the two blocks 
in the cold war, more evident in 1960. In general terms, a 
recurrent motif present in all years is the star-like struc- 
ture, formed by a central powerful economy surrounded 
by small dependent economies. The USA appears as one 
of those powerful hubs during all this period. However, 
other countries has gradually lost this role in favor of 
others. This is the case of the UK, which was the only 
star-like counterpart of the USA in 1948; since then its 
direct area of influence has been narrowing. On the con- 
trary, other countries have arisen for different reasons as 
new hub economies. This is the case of some European 
countries, Japan, and most recently, China. 



III. DIFFUSION ON COMPLEX NETWORKS 
AND THE DOLLAR EXPERIMENT 

The picture emerging from our analysis has intriguing 
similarities with other directed flow networks, such as 
metabolic networks [21j . that transport information, en- 
ergy or matter. Indeed, the trade imbalances network can 
be seen as a directed flow network that transport money 
In other words, we can imagine that net consumer coun- 
tries are injecting money in the system. Money flows 



along the edges of the network to finally reach producer 
countries. Producer countries, however, do not absorb 
completely the incoming flux, redistributing part of it 
through the outgoing links. The network is therefore 
characterizing a complex dynamical process in which the 
net balance of incoming and outgoing money is the out- 
come of a global diffusion process. The realization of 
such a non-local dynamics in the flow of money due to 
the trade imbalances spurs the issue of what impact this 
feature might have on the effect that one economy can 
have on another. In order to tackle this issue we per- 
form a simple numerical study, defined as the "dollar ex- 
periment" . The "experiment" considers running on the 
networks two symmetric random walk processes. Since 
we are limited by the yearly frequency of the empirical 
data, we assume at first approximation that the time 
scale of the changes in the structure of the underlying 
trade imbalances network is bigger than the characteris- 
tic diffusion time of the random walk processes. In the 
first case we imagine that a consumer country (As < 0) 
is injecting one dollar from its net debit into the system. 
The dollar at this point travels through the network fol- 
lowing fluxes chosen with a probability proportional to 
their intensity, and has as well a certain probability of 
being trapped in producer countries (As > 0) with a 
probability P a f, s — f^- . More precisely, if we consider a 
consumer country, such as the USA, the traveling dollar 
goes from country to country always following outgoing 
fluxes chosen with a probability proportional to their in- 
tensity. If in its way it finds another source it just crosses 
it, whereas if it finds a producer country j it has a prob- 
ability Pabs(j) of being absorbed. Mathematically, this 
process is a random walk on a directed network with 
heterogeneous diffusion probability and in the presence 
of sinks. By repeating this process many times it is possi- 
ble to obtain the probability that the traveling dollar 
originated in the source i is finally absorbed in the sink 
j. In other words, for each dollar that a source country i 
adds to the system, e%j represents the fraction of that dol- 
lar that is retained in country j. The symmetric process 
considers that each producer country is receiving a dol- 
lar and the traveler dollar goes from country to country 
always following incoming links backward chosen with a 
probability proportional to their intensity. If in its way 
it finds another sink it just crosses it, whereas if it finds 

a source j it has a probability P a bs{j) — of remain- 
ing in that country. The iteration of this process gives 
the probability gij that yields the fraction originated in 
the source country j of each dollar that a sink country 
retains. Consequently, these two quantities are related 
by the detailed balance condition 

\Asi\eij = Asjgji. (6) 

The matrices and g^ are normalized probability dis- 
tributions and, therefore, they satisfy that Y^j-sink e «i 
1 and J2i-,source9ji = 1- Using this property in the de- 
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TABLE II: Rankings from the Dollar experiment. Top: effect of two major source countries, USA and Switzerland, on the 
rest of the world. The first list is a top ten ranking of countries according to ejj , where the index i stands for the analyzed 
source. The second list is the top ten ranking of direct bilateral trade measured as the percentage of flux from the source 
country, that is, e l ij Cal = Fij/s° ut . Bottom: major contributors to two major sink countries, Japan and Russia. The first list is 
a top ten ranking of countries according to g^ , i standing for the analyzed sink. The second list is the top ten ranking due to 
direct trade. In this case, the direct contribution is g\° cal = Fji/s\ n . Countries in boldface have no direct connection with the 
analyzed country. The values for and gij are obtained from the simulation of the dollar experiment described in the text 
using 10 6 different realizations for each country, for the year 2000. 
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tailed balance condition, we can write 
Asj = ^2 e ij\As l \ and \Asi\ = ^ 9ji&Sj- 

i:source j:sink 

Then, the total trade imbalance of a sink or source coun- 
try can be written as a linear combination of the trade 
imbalances of the rest of the source or sink countries, 
respectively. Therefore, by measuring ey, it is possible 
to discriminate the effect that one economy has on an- 
other or, with gij , to find out which consumer country is 
contributing the most to a producer one, in both cases 
taking into account the whole topology of the network 
and the inhomogeneities of the fluxes. The advantage of 
this approach lies on its simplicity and the lack of tun- 
able parameters. Indeed, all the information is contained 
in the network itself, without assuming any kind of mod- 



eling on the influences among countries. 

By using this experiment it is possible to evaluate for 
a consumer country where the money spent is finally go- 
ing. For each dollar spent we know which percentage is 
going to any other producer country and we can rank 
those accordingly. It is important to remark that in this 
case countries might not be directly connected since the 
money flows along all possible paths, sometimes through 
intermediate countries. This kind of ranking is there- 
fore different from the customarily considered list of the 
first neighbors ranked by magnitude of fluxes. The anal- 
ysis indeed shows unexpected results and, as it has been 
already pointed out in other works [5] applying other 
methodologies, a country can have a large impact on 
other countries despite being a minor or undirect trad- 
ing partner, see Table 2. Similarly, producer countries 
may have a share of the expenditure of non directly con- 
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nected countries resulting in a very different ranking of 
their creditors. As an example, for each net dollar that 
the USA inject into the system, only 9.3% is retained in 
China although the direct connection imbalance between 
these countries is 16.7%. Very interestingly, we find that 
Switzerland spend a large share of his trade imbalance 
in countries which do not have appreciable trade with it 
and are therefore not directly connected such as Japan, 
Indonesia, and Malaysia. The Swiss dollars go to these 
countries after a long path of trade exchanges mediated 
by other countries. By focusing on producer countries 
we find other striking evidence. While the first importer 
from Russia by looking locally at the ranking of trade im- 
balances is Germany, the global analysis shows that one 
third of all the money Russia gains from trade is coming 
directly or undirectly from the USA. In Table 2, we re- 
port other interesting anomalies detected by the global 
analysis. 

IV. CONCLUSIONS 

In summary, we have introduced a novel quantitative 
approach applicable to any dense weighted complex net- 
work which filters out the dominant backbones while pre- 
serving most of the nodes in the original connected com- 



ponent. We have also discussed the behavior of a coupled 
dynamical process, the dollar experiment, which unveils 
the global properties of economic and trade partnerships. 
In a globalized economy, we face ever increasing problems 
in disentangling the complex set of relations and causality 
that might lead to crisis or increased stability. Focusing 
on just the bilateral relations among country economies 
is a reductionist approach that cannot work in a highly 
interconnected complex systems. We have proposed the 
use of the trade network representation and mathemati- 
cal tools that allow to uncover some basic ordering emerg- 
ing from the global behavior and the inclusion of non- 
local effects in the analysis of trade interdependencies. 
Future work on this grounds might help in the assess- 
ment of world trade relations and the understanding of 
the global dynamics underlying major economic crises. 
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